Shedding light on non-equilibrium dynamics of a spin coupled to fermionic reservoir 
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A single confined spin interacting with a solid- 
state environment has emerged as one of the 
fundamental paradigms of mesoscopic physics. 
In contrast to standard quantum optical sys- 
tems, decoherence that stems from these in- 
teractions can in general not be treated using 
the Born-Markov approximation at low temper- 
atures. Here we study the non-equilibrium dy- 
namics of a single-spin in a semiconductor quan- 
tum dot adjacent to a fermionic reservoir and 
show how the dynamics can be revealed in detail 
in an optical absorption experiment. We show 
that the highly asymmetrical optical absorption 
lineshape of the resulting Kondo exciton consists 
of three distinct frequency domains, correspond- 
ing to short, intermediate and long times after 
the initial excitation, which are in turn described 
by the three fixed points of the single-impurity 
Anderson Hamiltonian. The zero-temperature 
power-law singularity dominating the lineshape 
is linked to dynamically generated Kondo corre- 
lations in the photo-excited state. We show that 
this power-law singularity is tunable with gate 
voltage and magnetic field, and universal. 

Even though quantum dots (QD) are commonly re- 
ferred to as artificial atoms, their physical properties can 
be substantially different from that of real atoms. In 
particular, a single electron spin confined in a QD is sub- 
ject to hyperfine interactions with the QD nuclear spin 
ensemble and exchange interactions with the Fermi gas 
that controls the charging state of the QD. For many 
applications, such as those aimed at using QD spins to 
represent quantum information, the nuclear spin ensem- 
ble and the Fermi gas of electrons are treated as reser- 
voirs that induce decoherence of the confined electron 
spin. Tremendous progress in understanding and con- 
trolling spin decoherence induced in particular by QD 
hyperfine interactions has been achieved In contrast, 
relatively few studies have addressed the real-time dy- 
namics of an optically excited QD electron coupled to a 
nearby fermionic reservoir (FR). 
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The role of coherent tunnel coupling between a QD and 
a FR has been investigated since the late 1990s in the con- 
text of low-temperature transport spectroscopy0,[l,|3,[l| 
and is known to lead to one of the most spectacular 
phenomena of many-body physics - the Kondo effect 0]. 
The conductance through a QD in the Coulomb block- 
aded regime is proportional to the local density of states 
(LDOS) of the QD level; the Kondo effect introduces a 
quasiparticle peak in the LDOS, the Kondo resonance, 
at the Fermi energy ep, which leads to enhanced linear 
conductance through the QD. In contrast, optical sig- 
natures of the Kondo effect in quantum dots have so 
far not been observed experimentally 0, Q- Prior the- 
oretical work on this topic is also scarce, yet point to 
novel signatures: Govorov and co-workers predicted^], 
using a simplified approach based on variational wave 
functions [l0(, that doubly charged QDs could exhibit 
optical resonances whose width depends on an energy 
scale that is determined by the exchange interactions. 
Helmes et al. used Wilson's numerical renormalization 
group (NRG) method to calculate the absorption and 
emission lineshapes of excitons in the presence of Kondo 
correlations (henceforth called Kondo excitons), finding 
power-law divergences near threshold [ill]. 

Here we analyze the non-equilibrium dynamics follow- 
ing a quantum quench of a QD spin that is coupled to 
a FR. In a quantum quench, parameters of the Hamilto- 
nian of an interacting quantum system are changed over 
a very short time scale. We show that an optical absorp- 
tion experiment naturally implements such a quantum 
quench, and the measured optical absorption lineshape 
reveals, in a uniquely direct way, how for intermediate 
and long timescales standard quantum optical techniques 
such as Markov and Born approximations fail to describe 
the dynamics. Since the absorption of a photon suddenly 
switches on an attractive potential that (under suitable 
conditions) favors Kondo correlations and spin screening, 
the dynamic emergence of such correlations with time can 
be probed in ways not possible in transport experiments. 
For short times (large detunings), charge fluctuations 
dominate and coupling to the FR can be described using 
the Born-Markov approximation 12] . For intermediate 
times (intermediate detunings) spin fluctuations domi- 
nate; since the spin excitations of the FR are long-lived, 
the dynamics in this regime is non-Mar kovian. For long 
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FIG. 1: Schematics of the absorption configura- 
tion. The fermionic reservoir in the back-contact (of e.g. 
a Schottky-diode structure) couples via charge tunneling to 
the localized electronic degrees of freedom of the quantum dot 
(see Appendix 1). Starting from a an uncharged QD at t — 0, 
optical excitation to a b neutral exciton state (X°) at t = + 
featuring a single electron in the e-level of the QD (along with 
a hole in the h-level) which may develop Kondo correlations 
with the fermionic reservoir in the long-time limit. 



times (small detunings), the orthogonality catastrophe 
physics dominates and leads to a complete breakdown 
of the Born-Markov approximation. Remarkably, the re- 
sulting highly asymmetrical and singular absorption line- 
shape directly maps out these three regimes, which we 
link to the three fixed points of the single-impurity An- 
derson Hamiltonian. We find that the power-law singu- 
larity in the lineshape is tunable with gate voltage and 
magnetic field, and universal. 

Model for a tunable, optically active quantum 
impurity the excitonic Anderson Model. We 
consider a QD, tunnel- coupled to a FR (see Fig. 1), 
whose charge state is controllable via an external gate 
voltage V g applied between a top Schottky gate and 
the FR (see Appendix 1 for details). Assume V s to be 
tuned such that optical absorption occurs via the so- 
called neutral exciton transition[l4] (AT ), in which an 
electron-hole pair is created in the localized s-orbitals of 
the QD's conduction- and valence bands (to be called 
e- and h-levels, respectively). For circularly polarized 
light with frequency wl propagating along the 2-axis 
of the heterostructure . the QD-light interaction is 
described by Hi, tx (e^/i^.e _lWLt + h.c), where e\/h\ 
create an electron/hole in the e- and h-levels respec- 
tively and a = —a G {+,—}. We model the system 
before/after absorption by the initial/final Hamiltonian 
ffVf = + Hc + H t , where 
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describes the QD, with Coulomb cost U for double occu- 
pancy of the e-level, n ca — eje^, and hole energy Eha (> 
0, on the order of the band gap). The e- level's initial and 
final energies before and after absorption, e^ a (a = i, f), 
differ by the Coulomb attraction [/ c h(> 0) between the 
newly created electron-hole pair, which pulls the final e- 



level downward, e% a — e ca — S a i U e h (Fig. lb). This stabi- 
lizes the excited electron against decay into the FR, pro- 
vided that e l ca lies below the FR's Fermi energy ep = 0. 

H c = J2ka £ k<7 c la c ka represents a noninteracting con- 
duction band (the FR) with half- width D = 1/ (2p) and 
constant density of states p (st) = pO(D — \ek\) per spin, 
while H t = ^/T/TipY JIJ {elc a + h.c), with c a = c fc<r , 
describes its tunnel-coupling to the e-level, giving it a 
width T. A magnetic field B along the growth-direction 
of the heterostructure (Faraday configuration) causes a 
Zeeman splitting, e ca = e c + ^<rg c B, e\ lrT — + ^ag^B 
(see Appendix 1; the Zeeman splitting of FR states can 
be neglected for our purposes, see Appendix 8.) We set 
/ib = h = fcg = 1, give energies in units of D = 1 
throughout, and assume T,Bcrc(7, U c h <C D <C e^a- 
A realistic set of parameters would be, e.g. T k, 1-10 meV, 
U ~ U c h ~ 15-25 meV (in general U D h is slightly larger 
than U, resulting in a "trionic redshift" of the order of 
w 5 meV), D « 30 meV, e^s ~ 1-3 eV, g c k -0.6-0.7, 
<?h « 1.1-1.2. 

We focus on the case where the e-level is essentially 
empty in the initial state and singly-occupied in the 
ground state of the final Hamiltonian, n l e ~ and n l c ~ 1. 
(Here n° = (n e ) a is the thermal average with respect 
to H a of n c = J2a n ccr-) This requires e\ a ^> T, and 
— U + T < £g CT < — F. The initial ground state, needed 
below, will thus be approximated by the free Fermi 
sea, |G)i ~ rie fc ^<£ F c kcr\V' dC ) > neglecting terms of or- 
der r/£g CT . In particular, some (but not all) parts of 
the text will focus on the case that H f represents the 
symmetric excitonic Anderson model (if f =SEAM), with 
e l crT — —U/2, for which n l c = 1 exactly. 

Time evolution of the charge and spin after a 
quantum quench induced by absorption. To gain 
intuition for how the system would respond to the sud- 
den creation of an e\h\ exciton at time t = 0, it is 
instructive to calculate the subsequent time evolution 
of the average charge n c (t) = (n c+ + n e —)(t) and spin 
rh e (t) = \{n e + — n e -)(t) of the e-level, where n CIJ i{t) = 

Tr(e- iHtt p f p e iHtt n ea ,) and p { p = etp i e CT /(l - n 1 ^) is 
a projected version of the initial density matrix, nor- 
malized such that n ecr (0) = 1. Fig. 2 shows a typi- 
cal result for T = and i? f =SEAM, obtained using 
time-dependent NRG[il| (see Methods and Appendix 3). 
The non-equilibrium dynamics following such a quan- 
tum quench shows two distinct time scales: (i) Fluctu- 
ations in both charge and spin set in around the time 
scale t ~ l/l^eo-l associated with virtual transitions of 
electrons between e-level and FR. Whereas the charge 
equilibrates (towards 1) shortly thereafter, (ii) the spin 
decays (towards ~ 0) much more slowly, on the scale 
t ~ 1/Tk, where T K = JWj2e^ e ^+ u ^/ ( - 2U ^ is the 
Kondo temperature [13, EU associated with H l . (For fi- 
nite temperatures, the time scale on which fh c (t) decays 
is min{l/ 7 Kor, 1/T K }, where 7Ko r = T/\n 2 (T/T K ) is the 
Korringa relaxation rate[5j .) The decay is due to spin- flip 
processes, mediated by electrons of opposite spin hop- 



3 




FIG. 2: Non-equilibrium time evolution of spin and 
charge of the photo-excited electron. Plotted is the 
non-equilibrium time-evolution of charge and spin degrees of 
freedom of the photo-excited electron after the sudden cre- 
ation of an e\h_ exciton at time t = 0. While fluctuations 
in both e-level's total charge n e and spin set in around the 
time scale 1/le^l, the equilibration of the spin-<r populations 
n eCT (i) and screening of the local spin rh e (t) sets in on the time 
scale 1/Tk- The deviations (by about 3%) of n e (oo), n eo -(oo) 
and m e (oo) from their expected equilibrium values (1, ^ an d 
0, respectively, for the case # f =SEAM depicted here), are 
known artifacts of time-dependent NRG [3, presumably due 
to the NRG discretization scheme, which is inevitably coarse 
at large energies (see Appendix 3). Here, t is measured in 
units of 1/D. NRG parameters: U = 0.1D, s l e = 0.75(7, 
el = -0.5(7, T = 0.03(7, T K = 5.9 ■ 10~ 6 T, T = 0, B = 0, 
A = 1.8, Kept states: 1024, a = 0.4. 



ping between e-level and FR, leading to non-Markovian 
dynamics because the bath remembers its spin state be- 
tween two e-level spin-flips. As a result a FR screening 
cloud builds up over time, which ultimately screens the 
localized spin into a singlet. Note that this in particular 
means that for t > 1/Tk the "reservoir" (FR) is substan- 
tially modified, implying the complete breakdown of the 
Born-Markov approximation. We will come back to this 
point further below. 

The time-evolution depicted in Fig. 2 could in principle 
be observed by a 7r-pulse excitation of the QD followed 
by polarization resolved detection of the photolumines- 
cence. However, the fingerprints of the non-equilibrium 
dynamics can be more clearly discerned by measuring the 
absorption lineshape of a continuous-wave laser field, as 
we show next. 

Absorption lineshape of a Kondo exciton. Ab- 
sorption sets in once wl exceeds a threshold frequency 
u>th = Eq — Eq, which is on the order of e l cr7 + e\ lS - (mi- 



nus corrections due to tunneling and correlations). By 
Fermi's golden rule the absorption lineshape at temper- 
ature T and detuning v = ojl — u>th is proportional to 

A a (y) = 2nJ2pln\f(n\ein i \ 2 S(cJ L ~El l + E i m ). (2) 

Here \m) a and E^ are the exact eigenstates and -energies 
of H a and p x m — e~ Em ' T /Z 1 the initial Boltzmann 
weights. For future reference, we note that Eq. ([2]) can 
be expressed as A a (v) — — 2~lxaS° e {v), where, for T = 0, 

e c CT »=i(G|e CT ^4|G>i, (3) 

with v + = v + iO and H ! = H f — Eq — oj t h- Moreover, the 
Fourier representation G^{u) = J dte^^+^G^t), 

where G* (t) = -i6>(i)i(G|e iff Ve- lfftt 4|G)i, makes ex- 
plicit that this correlator directly probes the dynamics, 
described above and in Fig. 2, of a photo-generated elec- 
tron coupled to a FR (see Appendix for details). 

We used NRG to calculate A a {y) from Eq. @, gener- 
alizing the approach of Ref. [ll[ to T ^ by following 
Ref. [19| (see Methods and Appendix 2 for details). For 
clarity, we focus first on ff f =SEAM with 5 = 0. Fig. 3a 
shows a typical result: As temperature is gradually re- 
duced, an otherwise symmetric lineshape develops into a 
highly asymmetric one, dramatically increasing in peak- 
height as T — > 0. At T = 0, the lineshape displays a 
threshold behavior, vanishing for v < and diverging as 
v tends to from above. Fig. 3b analyzes this divergence 
on a log-log plot, for the case that T, which cuts off the 
divergence, is smaller than all other relevant energy scales 
(hence T = in all analytical calculations below) . Three 
distinct functional forms are discernible in the regimes of 
"large", "intermediate" or "small" detuning, labeled (for 
reasons discussed below) FO, LM and SC, respectively: 

(FO) | £ f ea | < v < D : A oc i/- 2 0(i/ - |4 CT |) ; (4a) 
(LM) T K < v < |4J : 4«-' \n- 2 {u/T K ); (4b) 
(SC) T < v < T K : A oc v~^' . (4c) 

A central goal of this paper is to explain the remarkable 
series of cross-overs described above. To this end we note 
that absorption at large, intermediate or small detuning 
probes excitations at successively smaller energy scales, 
corresponding to ever longer time scales after absorption, 
for which H can be represented by expansions H * + H' r 
around the three well-known fixed points[2(| of the AM: 
the free orbital, local moment and strong- coupling fixed 
points (r = FO, LM, SC), characterized by charge fluctu- 
ations, spin fluctuations and spin screening, respectively. 

Large and intermediate detuning dependence of 
the lineshape perturbative regime. For large de- 
tuning, probing the time interval t < 1/|4 CT | immediately 
after absorption, the e-level appears as a free, filled or- 
bital perturbed by charge fluctuations, described by the 
fixed point Hamiltonian H^ = H c + Hq D and the rel- 
evant perturbation H^ = H t . Intermediate detuning 
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FIG. 3: NRG vs. analytic results for the B — absorption lineshape. Solid lines represent the absorption lineshapes for 
_ff f =SEAM calculated by NRG a for several temperatures on a semi-log plot (See Appendix 5), and b at very low temperature 
(blue) on a log- log plot (inset shows the same lineshape in a linear plot). The latter reveals three distinct functional forms (red 
dashed lines) for high, intermediate and small detuning, labeled FO, LM and SC, respectively as captured by the fixed point 
perturbation theory (Eqs. ([5]), © and Q) (these, and thin blue line, were calculated assuming F 1 = 0). Arrows and light 
yellow lines indicate the crossover scales T, T K and |eL|- NRG parameters: U = 0.1D, 4 = 0.75(7, e* = -0.5(7, T = 0.03C7, 
T K = 5.9 ■ 10" 6 r, T = 3.3 ■ 10~ 10 r for b, B = 0, A = 1.8, Kept states: 1024, a = 0.4 for a and a = 0.5 for b. 



probes the times 1 / | < t < 1/Tk for which real charge 
fluctuations have frozen out, resulting in a stable local 
moment (assuming n ! c ~ 1); however, virtual charge fluc- 
tuations still cause the local moment to undergo spin fluc- 
tuations, which are not yet screened. This is described 
by [13, Hil Hl M = H c + const, and the RG-relevant per- 
turbation i?L M = j^s c ■ s c (a potential scattering term 
in i?L M , being RG-irrelevant, will be neglected in the 
discussion of the intermediate-detuning regime). Here 

= | Xw ijw V ( for 3 = e ' c )' are spin-operators 
for the e-level and conduction band, respectively (r are 
Pauli matrices), and J = 2UT / 1 7rg£ (e* +U) \ is an effective 
dimensionless exchange constant 18] . Constant contribu- 
tions to H* will not be specified, since they affect only 
Loth, whose precise value is not of present interest. It suf- 
fices to note that for both r = FO and LM, eJ.|G)i is an 
eigenstate of H* with eigenvalue w t h + Eq (within the 
accuracy of H*). 

For r = FO and LM, A a {v) can be calculated us- 
ing fixed point perturbation theory in H' r) as outlined 
in Methods and Appendix 5. At T = 0, we find 

4F 

A l °M = (\el\<v<D), (5) 

A™{») = f ^ (Tk<^<|4.I). (6) 

The crossover from LM to FO shows up in our numerical 
results as a small shoulder or side peak in A a (y) at | , 
depending on the detailed choice of parameters (easily 



seen from A^ M (|4 CT |) ^ ^°(kLl))- Tnis cross-over re- 
flects the fact that for v > \e f ea \, transitions from the h- 
level into unfilled states of the FR are possible, mediated 
by -ffpQ and using the e-level (Lorentzian-broadened by 
charge fluctuations, hence the v~ 2 dependence) as inter- 
mediate state, without creating additional particle-hole 
excitations. For intermediate detuning (LM) this is not 
possible; instead, absorption into the e-level is accompa- 
nied by the f/^ M -induced creation of electron-hole pairs 
in the FR, yielding an additional phase space factor oc v. 
In Eq. ([6]), we evoked scaling arguments [H[ to replace 
the bare exchange constant J by its renormalized, scale- 
dependent version, J(i/) = 1h _1 (i//Tk). Eqs. © and © 
reproduce Eqs. (|4"a|) and (|4b|). and quantitatively agree 
with the NRG results of Fig. 3. 

For clarity, the above discussion was confined to 
i/'=SEAM. However, it can be generalized straightfor- 
wardly to the non-symmetric case with e\, a ^ — |f7, as 
long as H l remains in the LM-regime, with fi l c ~ 1 (see 
Appendix 4). In the LM regime, the lineshape depends 
on e e(T and U only through their influence on Tk, and 
hence A^ M (v) is a universal function of v and 7k. This is 
illustrated in Fig. 4 for five lineshapes, shown in the lower 
left panel, corresponding to different choices of e ca and 
hence different TK-values (as indicated in inset). When 
these lineshapes are rescaled as A a {v) / '^(Tk) vs. u/Tk 
(main panel) , they collapse onto a universal scaling curve 
within the LM regime Tk ^ v < \£ l ca \- An experimental 
observation of such a scaling collapse would be a smoking 
gun for the existence of Kondo correlations. 
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FIG. 4: Universality in the LM-regime. Lower left panel: 
Five lineshapes, corresponding to five different choices of e eCT , 
indicated by color-coded dashed lines in upper inset, and ar- 
rows in main panel. The inset also gives fi e , n C e , An c and Tk 
as function of e ea . Main panel: When appropriately rescaled, 
the lineshapes collapse onto a universal curve in the LM- 
regime T K < v < |4 CT |- In the SC-regime T < v < T K , 
the curves do not collapse, since their exponents r\ a depend, 
via An c , on s etr . NRG parameters: U = 0.1D, U e h = 1.25(7, 
T = 0.062(7, T = 1.6 • 10~ 5 r, B = 0, A = 1.8, Kept states: 
1024, a = 0.5. 



Small detuning dependence and the tunable 
Kondo-edge singularity the non-perturbative 
regime. As v is lowered through the bottom of the LM 
regime, 3(v) increases through unity into the strong cou- 
pling regime, and A a (y) monotonically crosses over to SC 
behavior. In this limit, the assumption of a FR that re- 
mains unaltered by the interactions with the QD breaks 
down, requiring an alternative approach for determining 

For small detuning, i.e. long times t > l/"7k, a screen- 
ing cloud builds up that tends to screen the local mo- 
ment into a spin singlet, as visible in the decay of 
fh e (t — > oo ) — > (see Fig. 2). The screened spin singlet 
acts as a source of strong potential scattering for other 
FR electrons, causing the phase of each mode ka to shift 
by S a (£k<j) relative to its value for H 1 . This regime can 
be described by a strong-coupling fixed-point Hamilto- 
nian -ffgc + -^sc due ^° Nozieres (given in Appendix 6). 
It is formulated purely in terms of these phase-shifted 
c-electrons and makes no reference to e-level operators 
at all, since in the SC-regime the local moment is fully 
screened. Thus, the fixed-point perturbation strategy 
used above cannot be applied here. 

This hurdle can be overcome by working in the time- 
domain, and relating the correlator G° c (t) mentioned af- 
ter Eq. to the X-ray edge problem. The latter deals 



with the absorption lineshape for an incident X-ray to 
suddenly dislodge an electron from an atomic core level, 
placing it in the conduction band and leaving behind a 
core hole (i.e. a local scattering potential). Our situa- 
tion is similar, in that the screened singlet also serves 
as local scattering potential, but more complex, in that 
this potential is not turned on suddenly, but emerges 
only as the long-time limit of the dynamical build-up of 
a Kondo cloud to screen the local spin. Nevertheless, in 
the long-time limit an equation of motion approach can 
be used to relate the correlator G^ c (t) to a similar one, 
G£ c (t), which involves only conduction band electrons 
and whose form is known from the X-ray edge problem 
(see Methods and Appendix 6 for details). This readily 
leads to a power-law divergence characteristic of X-ray 
edge problems[22l H, 13, 

4°MaT K - 1 ( I //r K )-'' 1 (7) 

where the infrared singularity exponent [3, 

n a = 2An ra -^(A? w ) 2 , (8) 

a 1 

depends on the e-level's change in average occupation, 
An ca = n l ccr — figg.. Eq. ([5]) may be regarded as a gen- 
eralized version of "Hopfield's rule of thumb" [3, which 
was established rigorously in Refs. [H, [3- It nas an 
instructive physical interpretation, based on rewriting it 
as rja = 1 - E CT /(A< ff ,) 2 , where An' ca , = An 0(T - <W 
is the charge difference in level eer' between the final 
ground state |oo) (at time t — > oo) and the state |0 + ) 
the system finds itself in just after photo-excitation of 
a spin-cr electron (at t = + ). The "1" in n a repre- 
sents a v~ x power law divergence: it may be thought of 
as arising from a detuned, virtual transition into a nar- 
row e-level situated at v = (giving a Lorentzian detun- 
ing factor l/^ 2 ), followed by the creation of particle- hole 
pairs (with phase space v) to carry off the excess en- 
ergy v, resulting in a lineshape scaling as v/v 2 = v~ x . 
The X)cr'(A^LiJ 2 contribution to r\ a reflects Anderson 
orthogonality [13| : since |oo) and |0 + ) have localized ecr'- 
charges that differ by An' ca .,, their Fermi reservoir elec- 
trons see different scattering potentials, implying [3 that 
their overlap scales with effective system size L ~ v^ 1 as 
(oo|0+) -i-^-'( A "-') 2 . 

For the X-ray edge problem, An c<T and hence r\ a are 
fixed by material parameters (the strength of the core- 
hole potential). In contrast, in the present case they 
can be tuned experimentally by sweeping e\ a and e ea 
with a gate voltage or magnetic field (see Figs. 5b and 
Figure S3). This tunability can be exploited to study 
universal aspects of Anderson orthogonality physics that 
had hitherto been inaccessible. In particular, if the sys- 
tem is tuned such that n l c — Q and n l c = 1, Eq. can 
be expressed as rj a = \ + 2m l c a — 2(m i e ) 2 , where the final 
magnetization m f c — ^(n f c+ — n l c _) is a universal func- 
tion of g c B/TK- (At very large fields, however, a bulk 
Zeeman field, neglected above, will spoil universality, see 
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FIG. 5: Magnetic-field dependence of lineshape. a, Depending on whether the electron excited into the lower or upper 
Zeeman-split e-level (a = lower/upper, solid/dashed lines) the near-threshold divergence, A a (i>) cx v~ v ° , is either strengthened, 
or suppressed via the appearance of a peak at v ~ |g e -B|, respectively, b, Universal dependence on g c B /Tk of the local moment 
mi, calculated using the same values as in a for U and ei (dash-dotted line), and the corresponding infrared exponents rji OWCI 
(solid) and Supper (dashed) predicted by Eq. JS). Symbols: ^-values extracted from the near-threshold v^^" divergence of 
A a {y), for several magnetic fields (same color code as in (a)) and three values of V (distinguished by symbols). Symbols and 
lines agree to within 1 %, confirming the applicability of Hopfield's rule. NRG parameters: U = 0.1D, ej, = 0.75C7, = —0.5(7, 
r = 0.062(7, T K = 3.7 ■ 10~ 3 r, T = 0, A = 2.3, Kept states: 1200, a = 0.6. 



Appendix 8.) Thus, the exponents rj^^B) are universal 
functions of g c B/T^, with simple limits for small and 
large fields (see Fig. 5b): 



Slower/upper 



3 (1^1 «?k), 
±1 (|B|»T K ). 



(9) 



The notation a = "lower" or "upper" distinguishes 
whether the spin of the photo-excited electron (selectable 
by choice of circular polarization of the incident light) 
matches the spin of the lower or upper of the Zeeman- 
split e-levels, respectively. The sign difference between 
flower and Supper for \B\ 3> Tk arises since the change 
in local charge becomes fully asymmetric, An e , lower — ► 1 
while Art CiUppcr — > 0; as a result, for a = lower Anderson 
orthogonality is completely absent (Au^, = 0), whereas 
for a = upper it is maximal (An' ca , = 1). It follows, re- 
markably, that a magnetic field tunes the strength of An- 
derson orthogonality, implying a dramatic cr-dependence 
of the evolution of the lineshape A a (v) oc v~ ria with in- 
creasing \B\ (Fig. 5a): For ^4i owcr (^), the near-threshold 
singularity becomes stronger, tending towards v . In 
contrast, for ^4upper(^) the singularity becomes weaker, 
and once rapper turns negative, changes to an increas- 
ingly strong power-law decay, tending toward v +l \ this 
is accompanied by the emergence of an absorption peak 
near v = |<7 e -B|, associated with a transition into the up- 
per Zeeman-split level, broadened by Korringa relaxation 
of its spin (see Appendix 7). The fact that for \B\ > Tk 
a change of spin orientation of the photo-excited electron 
will turn a near-threshold divergence in the lineshape to 
a suppression that constitutes the low-frequency side of 
a broadened peak, is one of the most striking predictions 
of our analysis. 



Magnetic field dependence-dependence of the 
absorption threshold. The shift of the absorption 
threshold frequency wth = Eq — w ith magnetic field 
can be written as u)th(B) — Wth(0) = ^ag^B + 5u° h (B). 
The first term reflects the Zeeman energy of the photo- 
excited hole (which has pseudo-spin 3/2), the second the 
-B-dependence of the ground-state energy of the electron 
system. The general T = relation g c m% = 8E G /dB 
implies that the differential threshold shift offers a di- 
rect way of experimentally measuring the local moment 
difference between the final and initial ground states: 
d(5u e th )/dB = g e [m{{B) - m\(B)}. Also, for % ~ and 
h { c ~ 1, where Eq. (9) applies, this quantity can be re- 
lated to the infrared singularity exponents, allowing for a 
consistency check. Moreover, the asymptotic behavior of 
m l c for small fields (m e = —g e Bxo, where xo = 1/4Tk is 
the linear static susceptibility) and large field (|?rig| = h) 
implies: 



. ... _ , -(g c B) 2 /8T K 
^ h - 1 -g c B/2 



(\B\ «T K ) , 
(T K « B « |4|) 



(10) 



We note that a formula that interpolates through 
both of these regimes is given by Suj° h = Tk — 

\\J (9cB) 2 + (2Tk 2 ) which, up to numerical prefactors, 
has a similar functional form as the expression for the 
ground state energy of the s-d model derived in Ref. [28| . 
The quadratic B-dependence of <W{r h for small fields of- 
fers a straightforward way to determine the Kondo tem- 
perature experimentally. The accessibility of the in- 
dependence of the ground state energy and e-level mag- 
netization via the absorption threshold is a remarkable 
advantage of the proposed optical probe of Kondo physics 
in this paper - these quantities are not accessible via 
transport measurements. 
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We have demonstrated that optical absorption in a sin- 
gle quantum dot can be used to implement a quantum 
quench - a sudden change in the Hamiltonian governing 
the dynamics of the many-body system. Given that the 
relevant quantum dot parameters are tunable via exter- 
nal electric and magnetic fields, this system constitutes 
one of the rare, if not unique, experimentally accessi- 
ble solid-state systems where a tunable quantum quench 
could be realized. Our work sets the stage for exploring 
numerous further interesting problems, such as (i) the ef- 
fect of an exchange interaction between e- and h- levels; 
(ii) the effect of nuclear spins on the electron spin dynam- 
ics; (iii) using a pump-probe protocol to study the non- 
equilibrium time evolution even more directly; (iv) study- 
ing the coherence of optical Raman transitions, where the 
virtually excited intermediate state is a Kondo correlated 
state of the QD electron and an adjacent fermionic reser- 
voir, and (v) attempting to exploit strong correlations to 
realize a "Fermionic quantum bus" between two distant 
QD spins. The presented results have an important bear- 
ing on the quantum control of a spin degree of freedom in 
nano-structures and on quantum optical techniques used 
to study reservoirs composed of fermionic constituents. 



transform of B(t), B(ui) = J dt e luJt B(t), can be expressed 
in Lehmann representation: 

= f(n'|/5» ff (n|B|n')f ■ 2ir5 (w - E f n , + E f n )(ll) 



This expression can be calculated using FDM-NRG. See 
Appendix 3 for details. 



C. Fixed-point perturbation theory 

To calculate A a {v) at T = for v > T K , set H f -> 
H* + H' r in Eq. and expand in powers of H' r , for r = 
FO or LM. At zeroth order this yields a peak j4 ct (^) ~ 
8{v), which is irrelevant for v > Tk- To lowest non 
vanishing order in H' r 



A r a (u) = - Aim Ug^K ~K4 |G)i 

u z L v + - Hi 



which we shall evaluate for T = B = 0, |er* CT | - , 
|G>i^n £fc „< eF 4jVac). We find 



(12) 
hU and 



Methods 



Kiy) 



(13) 



A. Numerical Renormalization Group for 
calculation of optical absorption lineshape 

The optical absorption lineshape given by Fermi's 
golden rule Eq. |(3J), can be calculated at finite tempera- 
tures using full density matrix (FDM) numerical numer- 
ical renormalization group (NRG) [191 . Because Eq. © 
contains matrix elements between initial and final eigen- 
states of different Hamiltonians, H 1 and H l , two separate 
NRG runs (NRG run #1 and #2) are required to calcu- 
late the initial and final eigenstates ({|m)i} and {|n)f}) 
as well as eigenenergies {{E^} and {E^}). The double 
sum in Eq. |(5J), over all initial and final eigenstates, is 
performed via a "backwards" run from the end to the 
beginning of the Wilson chain (l6j: for each shell k, the 
contribution towards the initial density matrix p 1 from 
that shell (obtained using data from NRG run #1), and 
the matrix elements |f (n|ej.|m); | between shell-A; eigen- 
states from NRG runs #2 and NRG #1 are calculated, 
and binned according to the corresponding frequency dif- 
ference E f n — E x m . See Appendix 2 for further details. 



B. Non-equilibrium dynamics via NRG 

The expectation value of an observable B after ab- 
sorption is given by Bit) = Tr ^/5 p (i)_B^ where the time 
evolution is governed by the final Hamiltonian, p p (t) = 
e~ lH f * p t,e* H f * . For reasons discussed in Appendix 3 we 



where T™(u) = 2(T/irp)F v {-i6{t) i (G|c T (t)4(0)|G) i } 
in the FO-regime {T v is the Fourier-transform operator 
T v {A} = J dte lt ^ +ultb) A(t)). In this regime, the ab- 
sorption process can be understood as a two-step process 
consisting of a virtual excitation of the QD resonance, fol- 
lowed by a tunneling event to a final free-electron state 
above the Fermi-level. Therefore, the absorption rate is 
proportional to the tunneling density of states (into the 
FR) at detuning v. On the other hand, in the LM regime, 

T™{v) = 3/8(J/p) 2 FA-i(>(t)i(G\s° c (t)s-*(0)\Gh} 
where we used (s z c (t)s z c (0)) = (s+(t)s~ (0))/2. This re- 
sult reiterates that the intermediate detuning probes 
spin fluctuations, as observed in the dynamics (Fig. 2). 
Evaluating these correlation functions yields Eq. (f5|) and 
Eq. ©. In Eq. ©, we have inserted the scale-dependent 
exchange interaction J{v) — In - (z//Tk) which results 
from the logarithmic enhancement of the exchange inter- 
action within perturbative RG. The finite-temperature 
form (T 3> Tk) of the lineshape is discussed in Ap- 
pendix 5 (see also Fig. 



D. Absorption in the strong-coupling regime and 
Fermi-edge physics 

Eq. §5§ can be expressed as the Fourier transform, 
Geei v ) = {Gl c {t)}, of a correlator 



Gl e (t) = -i6(t)(e„{t)et)i 



(14) 



find it convenient to take p f p = p p /[l 



n l ca \. The Fourier 



involving operators defined to have an anomalous time 
dependence, 0(t) — e lH t Oe~ lH 1 . This anomalous time 
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dependence, involving both H 1 and H { , reflects the fact 
that the creation of a hole during optical absorption 
abruptly lowers the e-level. We next relate Qc C { v ) to a 
similarly-defined correlator of FR electrons, 

Gutty) = ?v {-**(*) Mt)cLh} (is) 

using equations of motion in the asymptotic limit v — > 0: 



£fe CT )(^+ + A - £ k ,„) ' 



(16) 



where A = uj t h — £\ia ■ This asymptotic relation is estab- 
lished in Appendix 6 and implies 



irpT 



(17) 



where G°Jy) = £ fcfc ,S&,(i/). To calculate G c CT c (t) for 
t ;§> l/Tk, we may now replace H 1 — > _ff c and i? f — > 



G-W-^GIe^'c^e-^c^tiG). 



(18) 



This response function is similar to that calculated in 
the X-ray edge problem [22I l23l . l24l ] and its calculation is 
standard (e.g. [2J|) and yields (we show only the leading 
power law) 



(19) 



where 8 a = 8 a {Q) denotes the phase shifts at the Fermi 
energy. The phase shifts at the Fermi energy, in turn, 
are given by i5 CT = 7rAn C(T , according to the Friedel sum 
rule [29(| . valid for T = and for arbitrary values of B. n l c 
and rig. Collecting results, we find Eq. Q and Eq. {SJ. 
See Appendix 6 for details. 



Acknowledgments 

AI and HET acknowledge support from the Swiss NSF 
under Grant No. 200021-121757. HET acknowledges 
support from the Swiss NSF under Grant No. PP00P2- 
123519/1. BB acknowledges support from the Swiss NSF 
and NCCR Nanoscience (Basel). JvD acknowledges sup- 
port from the DFG (SFB631, SFB-TR12, De730/3-2, 
De730/4-l), the Cluster of Excellence Nanosy stems Ini- 
tiative Munich and in part the National Science Foun- 
dation under Grant No. NSF PHY05-51164. AI ac- 
knowledges support from an ERC Advanced Investigator 
Grant, and LG from NSF Grant No. DMR-0754613. 



APPENDIX 1: THE MODEL AND ITS 
EXPERIMENTAL REALIZATION 

We focus here on differential transmission spectroscopy 
of single semiconductor QDs in gated heterostructures 



with an adjacent tunnel- coupled electron reservoir. A 
possible realization is a Schottky diode structure 14 
where a voltage applied between a top Schottky con- 
tact and a two-dimensional electron gas (2DEG) in a 
modulation-doped GaAs layer underneath a layer of self- 
assembled InAs QDs, is used to adjust the relative energy 
of the QD electron with respect to the Fermi-energy of 
the 2DEG (Figure 151)). To realize the situation in Fig. 1 
we consider here a range of gate voltages for which the 
QD is uncharged before the optical excitation (See Fig. 4, 
upper inset). In contrast to most of the earlier quan- 
tum optics experiments where the separation between the 
electron gas and the QD layer was kept large to ensure 
weak spin-flip co-tunneling, we focus here on structures 
with a small barrier where tunnel coupling of a QD elec- 
tron is only a factor of 10 weaker than its single-electron 
charging energy. We estimate that with a 15 nm neutral 
GaAs barrier, the strength of the exchange interactions 
would be strong enough to yield Kondo temperatures ex- 
ceeding 100 mK. 

Spin-orbit interaction in InGaAs QDs leads to a split- 
ting of about 200 meV between the J = 1/2 and J = 3/2 
valence band states. In the limit of asymmetric quantum 
confinement typical for our QDs, the four J = 3/2 bands 
further split into two Kramers doublets. The states with 
J z = ±3/2 along the growth (= strong confinement) di- 
rection have the lowest zero-point energy due to their 
heavy mass for motion along the growth direction and 
hence are the ones relevant for lowest energy optical ex- 
citations. The confincment+strain induced splitting be- 
tween the J z = ±3/2 (heavy-hole) and J 2 = ±1/2 (light- 
hole) states is estimated to be around 20 meV. In the 
main text, the spin label of the valence band electrons 
(or holes, h a ) therefore refer to pseudo-spin J z = ±3/2, 
and the hole Zeeman energy has the form a^g^B. 

This situation leads to near perfect correlation between 
the circular polarization of the excitation light and the 
spin of the optically excited electron in InGaAs QDs (as 
well as any QD structure with possibly the exception of 
chemically synthesized nanocrystals) . The hole g- factor 
in these QDs is highly anisotropic: along the growth di- 
rection, the value is around 1.2. In plane g-factor ranges 
from - 0.5. 

For the calculations in the main body of the paper, we 
assume the electron-hole recombination rate (typically of 
order r c h ~ 1/xeV) to be negligible compared to all other 
scales, take the hole level to lie within the gap of the FR 
and neglect its coupling to the latter. The newly-created 
electron-hole pair will experience a mutual Coulomb at- 
traction J7 c h(> 0) and a much weaker exc hang e interac- 
tion J h (typically of order J h ~ 200 /zeV) [15| . 

While the interplay of the intra-dot electron-hole ex- 
change interaction and the Kondo correlations is interest- 
ing on its own right, the effects of electron-hole exchange 
can experimentally be avoided in two ways. One possi- 
ble alternative is to start out with a single-hole charged 
QD; the experiments on single QDs embedded in n-type 
Schottky structures have already shown that the life- 
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a b 




FIG. SI: The Schottky diode structure for the exper- 
imental observation of the optical signatures of the 
Kondo effect. 

time of such an optically charged state could well exceed 
100 fisec. Resonant optical excitation of such a single- 
hole charged QD leads to the formation of the so-called 
X 1+ trion which consists of two valence-band heavy- holes 
forming a singlet and a single e-level electron with van- 
ishing electron-hole exchange interaction. Another possi- 
ble alternative is using indirect excitons in coupled QDs 
where the electron and the hole wave-functions have van- 
ishing spatial overlap. 

Here, we assume J c h <C T%, or an experimental re- 
alization where it can neglected, so that the hole has 
no dynamics after creation. With the hole "frozen", 
the electron dynamics can be described by a pure An- 
derson model (AM), Eq. (1), for the e-level (involving 
only e CT , not h s ), whose initial and final energies be- 
fore and after absorption, e" ff (a — i,f), are related by 
e" CT = e ca — 5 a {U e h (Fig- lb). The downward pull of the 
Coulomb attraction U c h stabilizes the excited electron 
against decay into the FR, provided that e ea lies below 
the FR's Fermi energy e-p = 0. 



APPENDIX 2: NUMERICAL 
RENORMALIZATION GROUP AND FERMI'S 
GOLDEN RULE 

The numerical renormalization group [13] is an itera- 
tive method for numerically diagonalizing quantum im- 
purity models such as Anderson impurity Hamiltonians 
H a (a = i, f) specified around Eq. (1). The spectrum of 
states of the Fermi reservoir is coarse-grained using a log- 
arithmic discretization scheme governed by a parameter 
A > 1 (typically A — 2), followed by an exact mapping 
of the discretized model onto a semi-infinite chain, the 
so-called Wilson chain, whose hopping amplitudes decay 
exponentially along the chain, as ifc ~ A _fc / 2 . This pro- 
duces a separation of energy scales and makes it possible 
to diagonalize the Hamiltonian iteratively: knowing the 
eigenstates of a chain of length k — 1, one adds site k and 
calculates the "shell" of eigenenergies of the Hamiltonian 
for the chain of lenght k. The high- lying eigenstates of 
that shell are "discarded" , while the low-lying states are 
"kept" and used for the next iteration. The spectrum of 
eigenenergies so obtained typically flows past one or more 
non-stable fixed-points and finally converges towards a 



stable fixed point, whereupon the iterative procedure can 
be stopped. In practice one thus deals with a finite Wil- 
son chain, whose length is set by the smallest energy 
scale in the system (e.g. the Kondo temperature, tem- 
perature, or magnetic field). By combining NRG data 
from all iterations, it is possible to construct a complete 
set 16] of approximate many-body eigenstates of the full 
Hamiltonian. These can be used to evaluate equilibrium 
spectral functions via their Lehmann-representations; at 
finite temperatures, this can be done using the full den- 
sity matrix (FDM)-NRG[l|. 

Since Eq. (2) expresses the Fermi golden rule absorp- 
tion rate via a Lehmann representation, it, too, can be 
evaluated using NRG [Til]. However, it contains matrix 
elements between initial and final states that are eigen- 
states of different Hamiltonians, H 1 and H f . Hence, two 
separate NRG runs are required to calculate these (simi- 
lar in spirit to what is done for time-dependent NRG[16]). 
The strategy is then as follows: 

• NRG run #1 generates a complete set of approxi- 
mate eigenstates |m)i and eigenenergies E l m for the 
intial Hamiltonian H 1 (without exciton). 

• NRG run #2 generates a complete set of approxi- 
mate eigenstates \n)f and eigenenergies for the 
final Hamiltonian H (with exciton). 

• The double sum in Eq. (2), over all initial and final 
eigenstates, is performed via a "backwards" run, 
with site index k running from the end to the be- 
ginning of the Wilson chain [16]: for each shell k, 
the contribution towards the initial density matrix 
p 1 from that shell (obtained using data from NRG 

run #1), and the matrix elements f (n|ej.|m)i| be- 
tween shell-fc eigenstates from NRG runs #2 and 
NRG #1 are calculated, and binned (see below) ac- 
cording to the corresponding frequency difference 

E li ~ E ln- 

• The T = threshold frequency for the onset of ab- 
sorption is given by the difference of groundstate 
energies of NRG runs #2 and #1, w t h = Eq - E h- 
The absorption spectrum is expected to have di- 
vergences at the threshold w t h, hence all frequency 
data are shifted by the overall threshold energy uth 
prior to binning. (For finite temperature, the sharp 
onset is broadened and divergencies are cut off.) 

• The discrete eigenenergies of shell k are spread 
over an energy range comparable to the character- 
istic energy K~ k / 2 scale of that iteration, which de- 
creases exponentially with k. Thus, the bins used 
for collecting the discrete data are likewise chosen 
to have widths decreasing exponentially with de- 
creasing energy. The discrete, binned data are sub- 
sequently broadened using a log-Gaussian broad- 
ening scheme, characterized by a broadening pa- 
rameter a 19], typically taken as a = 0.6. For fi- 
nite temperature, the maximum of the absorption 
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peak occurs at a frequency ^ max , that is slightly 
larger than the T = threshold frequency at v = 
(see Fig. 3a). Thus, for finite temperature, the bin- 
ning procedures (both the binning frequencies and 
the frequency ranges in which we change from log- 
Gaussian to Gaussian) for frequencies above and 
below t'max were set up to be symmetric with re- 
spect to !/ mm . This ensures that the A a (y) curves 
for v smaller and larger than y max match smoothly. 

APPENDIX 3: TIME EVOLUTION VIA NRG 

The nonequilibrium time evolution of the system in re- 
sponse to the sudden creation of an e\h^ exciton at time 
t = can be calculated using time-dependent NRG [IB]. 
The sudden addition of a cr-electron in the e-level projects 
the system's initial density matrix p 1 onto a projected 
density matrix p p , 

p' = Yl \ m h(p x m) iH Pp = 4p' e <r ■ (si) 

m 

Its norm is smaller than 1, 

Trp p = Tr ) = 1 - < CT , (S2) 

reflecting the fact that states in which the e-level had 
already contained a cr-electron prior to absorption are 
projected to 0. (n 1 ^ = (h crT )i denotes an expectation 
value w.r.t. to H 1 .) After this projection, the e-level oc- 
cupancies for spin a' are reset (or "reinitialized" ) to be 

TV,* * ^ _ / 4,oo + 4,o<t = 1 - 4<r for g' — g /qo\ 
{Ppea) ~\nl m {oxa' = a^ 6) 

where n], 00 and n\ os are the initial probabilities (with 
respect to H 1 ) for the e-level to have been completely 
empty (00), or to have contained no cr-electron but an 
a-electron (Oct), respectively. 

For the sake of studying the time evolution after ab- 
sorption, we find it convenient to take p p = p p /[l — nJ, CT ] as 
starting density matrix, normalized such that Trp p — 1 
and Tr p p n ca = 1. Its subsequent time evolution is gov- 
erned by the final Hamiltonian, p p (t) = e~ lH t p t p e lH l . 
The corresponding expectation value of an observable B 
is thus given by 

B(t) = Tr (g{t)B) . (S4) 

Fourier-transformed to frequency space using B(u>) — 
J dt e iwt B(t), one obtains 

B(w) = t(n'\p p \n) f{ (n\B\n') t ■ 2nS ( u - E l n , + £*)(S5) 

n,n' 

The latter expression, being written in Lehmann repre- 
sentation, is again well-suited for FDM-NRG [H . Follow- 
ing [l6l ]. it is again possible to formulate the procedure 



such that matrix elements are always calculated within 
the same Wilson shell. 

The discrete data contributing to B(ui) is binned and 
smoothened using a log-Gaussian broadening function, 
as described in the previous section, but using a smaller 
broadening parameter, a — 0.3. B(t) is obtained by 
Fourier-transforming the broadened version of B(u>) to 
the time domain at the very end of the calculation. 

The values obtained for n e (oo), n ea (oo) and m e (oo) de- 
viate from their expected equilibrium values (1, i and 0, 
respectively, for the case 7J f =SEAM depicted in Fig. 2) 
by about 3 %. This is a known artifact of time-dependent 
NRG [l|| , presumably due to the fact that the NRG dis- 
cretization scheme inevitably is rather coarse at large en- 
ergies. 



APPENDIX 4: PERTURBATION THEORY FOR 
THE EXCITONIC ANDERSON MODEL 

The optical absorption lineshape for a g — + excitation 
and the excitonic Anderson Model can be calculated via 
a perturbation expansion up to second order in the tun- 
nel coupling v — y/T /irp (we keep here the momentum 
dependence of Vk) for the initial and final states entering 
Eq. (2). (The wavefunctions used below are similar in 
structure to those used in Ref.[l3] for variational calcu- 
lations of various emission and absorption spectra of Ce 
compounds.) 

We define the empty e-level, h-level and unperturbed 
FR as |0) = |0) o (g> |0)h ® \0)f- The bare excitonic reso- 
nance frequency is loq = Wth = £c — U c \l + £h which also 
marks the approximate onset of the edge (threshold) be- 
havior. We write the initial state |G)i expanded to second 
order in Vk as 



a 

+ E (F " k t Fl 4>^]|0) (S6) 

fc'>0 

The set of final states are defined by particle number 
conservation. A few possible zeroth order final states 
accessible from the initial state (to low orders in Vf.) are: 
c k a h^\0) (X + exciton), etcfc^ejftt |0) (X~ exciton) and 

% u^k' <?&}J%\fy + e_n excitation). Note that the 
last class of states includes for k = k' the bare neutral 
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exciton final state. The relevant final states are then 



\xt) 



1 



E 



£fc — £ c 
Vk' 



£fc' 



T Ck> 



eh 



0>) 



l**> 



£ c + U - U ch - E k 



E 



£ e + f/ - f/ eh - £fe' 



"fe,T e T 



(S7) 

ejcfc.iejHl ) 
(S8) 



l^fcerfc'er'/j) 



1 + 5* , 



Vk 



£k 



+ 6 a > 



C/ + C/ch 



,Ck,-v 



(£fc-£fc')(£c-^oh-£fc') ' ' 



^-mE 



(£fe - £g)(£fc - £c - £7 + U eh ) 



4<Av4Hl°> ( S9 ) 



We have included the lowest order contributions in c-h 
excitations in the FR. The final state energies to lowest 
order in Vk are given by 



E, 



(x+) _ 



Wth + [£fe - (£ c 



U, 



ch) 



Ef~ ] = wth + [(£c - C/ch + U)- £fe] 



fx ) 

^0 = ^th + £fe - £fe' 



(S10) 

(Sll) 
(S12) 



For SEAM, the analytical form of the lineshape is: 
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FIG. S2: The absorption lineshape at high tempera- 
tures. For T ^> 7k, the log-log plot shows the two portions 
of the absorption lineshape for red-detuning iv < /v max ) and 
blue-detuning (v > v mam ) with respect to the frequency /v max 
(at which A a (y) reaches its maximum). Here, 7Kor = 0.15 T. 
NRG parameters: U = 0.1D, ej, = 0.75(7, e{ = -0.5(7, 
T = 0.03(7, T K = 5.9 ■ 10~ 6 R T = 100 Tk, B = 0, A = 1.8, 
Kept states: 1024, a = 0.5. 



APPENDIX 5: FINITE TEMPERATURE 
FIXED-POINT PERTURBATION THEORY 

To calculate A a {v) for T > T K set H l -> H r * + i?; in 
Eq. (3) and expand in powers of H' r , for r = FO or LM. 
Note that the proper finite-T generalization of Eq. (3) is 
Qeei 11 ) = 3~v {~ i9(t) ( e <r(i)e^)i} . Let us focus only on the 
temperature dependence of the lineshape for \v\ < |e eo .| 
and T>T K : 



3tt/4 7Kor (^,r)/7r 



(S14) 



where 7koi(^, T) is the scale-dependent Korringa relax- 
ation rate given by 



7Kor(^,T) 



^r/ln 2 |T/T K | forH<T, , . 
7T7/ /ln 2 |iVT K | for >T 1 ' 



This expression is compared to the lineshape calculated 
by NRG in Fig. [51 Note that for v > T, we recover 
Eq. (6). 



The singularity at v = U/2 can be related to a 
tunneling- assisted X~~ resonance. The intermediate state 
there is an X" charging state which is reached from the 
X° state by a FR electron tunneling in, and thereby sat- 
isfying energy conservation. We expect that a calculation 
including higher orders of perturbation theory will regu- 
larize this singularity. 

The FO "threshold" term Q(y — -j), which corresponds 
to Eq. |5]). is specific to the symmetric Anderson model. 
For the asymmetric Anderson model, this threshold splits 
into two, at v = \e{„\ and v = el„ + U. 



APPENDIX 6: EVALUATION OF THE 
ABSORPTION LINESHAPE IN THE 
STRONG-COUPLING REGIME 

As we descend to detunings well below the Kondo 
scale, v <C Tk, the physics is governed by the strong 
coupling fixed point. It describes a fully screened sin- 
glet, acting as source of strong potential scattering for 
other FR electrons, causing the phase of each mode ka 
to shift by 8 a (eka) relative to its value for H 1 . According 
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to Nozieres[2l| (see also [5]), the fixed point Hamiltonian 



at B = is Hg C — J2ka £ kcl<jCka' wnere tildes denote 
operators representing phase-shifted modes. The lea ding 
relevant perturbation for < B < Tk has the form[5ll2lf 



a^ka 



■p a /{-KpT K ) 



We shall not explicitly use this fixed point Hamilto- 
nian, however, since the strategy (described in Methods 
section) of perturbing around the fixed point is of no use 
for calculating G° c {v) of Eq. (3). The reason is that G° e {v) 
is formulated in terms of e ca and eJ CT operators, whose 
dynamics is determined by higher-energy excitations of 
the FR not described by Hg C + H§ c . To circumvent 
this problem, we use an equation of motion approach to 
first derive an asymptotic relation between the impurity 
Green's function and the Green's function of the itinerant 
electrons. 

We start by noting that the correlator occuring in 
Eq. (3) can be expressed as the Fourier transform, 
Qeeiy) = {Gec(*)}j °f a correlator 



GIM = -i6(t) <e ff (t)4)i , 



(S16) 



involving operators defined to have an anomalous time 
dependence, 0(t) — e lH t Oe^ lH *. This anomalous time 
dependence, involving both H 1 and H { , reflects the fact 
that the creation of a hole during optical absorption 
abruptly lowers the e-level. 

To relate Ge e (v) to a similarly-defined correlator 
G% c {v), we n °t e that the anomalous time-dependence of 
Cka(t) implies the equation of motion 

idtc k „{t) = [ck*(t),H i }+ Ck*(t)(e hs - U ch n c ) , (S17) 

= Cko{t)[Ska + Ehff - U ch n c ] + Ve<r(t) , (S18) 



where v = y/T/np. Inserting this into the definition of 
G a kk ,{t) one finds 

id t G% k , (t) ~ [e ka + e h9 ]G% k , It) + vG° k , (t) , (S19) 

where terms that become subleading for t — > oo have 
been dropped (a term containing 5(t); and one containing 
{ckcr(t)n c c kl(T )i which contains more operators and hence 
decays more quickly with time than the correlators that 
were retained). Similarly one finds 

id t GZ k , (t) ~ [e w<r + e hs ]Gl k , (t) + vG° c (t) , (S20) 

where the cyclic property of the trace was used to write 

<[e eCT (i),ir]4 v >i = (c ca {t)[H\cl, a ]) { . (S21) 



Fourier-transforming Eqs. (|S20[) and (|S19[) using the con- 
vention stated just before (|S16[) and eliminating G° k ,(v), 
we readily find the asymptotic relation: 



v 2 GM 



, (S22) 



(v+ + A - e ka ){v + + A - sk'a) 
where A = u>th — It follows that 

= E ^ W ~ -«p t GM > ( S23 ) 



where the double sum Y^kk' = P 2 I d e fcdefc' is exhausted 
by two 5-functions, since |A| is of order |e e(7 |, which we 
assume to be smaller than the bandwidth D. Eq. (|S23|) 
implies 



irpT 



(S24) 



To calculate G° c (t) for t » l/?k, we may now replace 
H 1 -> H r and H f 



#sc + #sc in Eq- (H2l : 



(t) ~ i(Gje l 



-iHS c tJ 



|G)i 



(S25) 



This response function is similar to that calculated in 
the X-ray edge problem [12, HH [24[ : there absorption of 
an X-ray photon excites an atomic core electron into the 
conduction band (described by cj.), leaving behind a core 
hole which constitutes a scattering potential with respect 
to H c (described by H$ c ). The calculation of Eq. (|S25|) is 
standard (e.g. [HJ|) and yields (we show only the leading 
power law) 



G° c (t) ~ t -iv.-*)'+ny*> 



(S26) 



where S a = S a (0) denotes the phase shifts at the Fermi 
energy. This power-law has an instructive interpretation, 
due to Hopfield 2q| : according to Anderson's orthogo- 
nality catastrophe [13, two Fermi seas subject to differ- 
ent local scattering potentials that cause their modes to 
differ in phase by 5 CT (efe CT ), have a ground state overlap 
f{G\G)i ~ L^- 5 ^ )/^ 2 which vani shes in the limit of 
system size L — > 00. In analogy, Eq. (|S26[) can be viewed 
as the overlap f (G|cJ.|G)i for systems of size L oc t: the 
effect of cj., which puts an extra spin-cr electron at the 
scattering site at t = 0, is analogous to having an addi- 
tional infinitely strong scatterer of cr-electons in the ini- 
tial, but not the final state, implying an extra shift — ir 
for the phase difference of the S a modes. 

The phase shifts at the Fermi energy, needed in 
Eq. (|S26p . are given by 8 a = 7rAn C(T , according to the 
Friedel sum rule[Ii| [3C|, valid for T = and for arbi- 
trary values of B, n f c and rig. Collecting results, we find 
from Eqs. (jS24f and (|S26|) that 



4 C (-) 



(S27) 



with the infrared singularity exponent r\ a given by 
Eq. (8), 



77 CT = 2A7i CCT -E( A? w) 2 



(S28) 



The dimensionful prefactor in Eq. (|S27[) was adjusted 
to ensure that Eqs. ([7]) and ^ match, up to numerical 
prefactors, at the crossover scale v — Tk, implying a pref- 
actor (pT)T^^ 1 in Eq. (fS26|) . We have checked Eq. (8) 
numerically for a range of parameter combinations, see 
Fig. 5b and Fig. IS31 finding it to hold to within 1 %. 
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FIG. S3: Checking Eq. (8), Hopfield's rule of thumb 
for the singularity exponent r\ a of the Kondo-exciton, 
for various combinations of |e e | and B. Solid lines were 
calculated via the the r.h.s. of Eq. (8), using NRG-results for 
An ccr as input; symbols give the values of the exponent r/ CT 
extracted from A a (y — > 0). Symbols agree with solid lines to 
within 1 %, confirming the validity of Hopfield's rule of thumb 
for absorption into a Kondo-correlated final state. NRG pa- 
rameters: U = 0.1D, r = 0.1(7, T = 0, A = 2.3, Kept states: 
1700, a = 0.6. 



APPENDIX 7: MAGNETIC FIELD 
DEPENDENCE OF THE ABSORPTION 
LINESHAPE 

In the main text, we pointed out that as \B\ in- 
creases past Tk, the evolution of the infrared exponents 
r\ a with B depends on whether the spin a of the photoex- 
cited electron matches that of the lower or upper of the 
Zeeman-split e- levels, which we distinghuished by writ- 
ing a = lower or upper, respectively. Here we discuss this 
cr-dependence in some more detail. 

It is instructive to express r\„ of Eq. (8) in terms of 
the final occupation n l e and magnetic moment m f c of the 
e-level, writing n l crT = ^n l c + am : 

7 h = % + 2m{ a - 2(m f e f , m = n*(l - ifi*) . (S29) 

In particular, this implies the simple relation m e = 
j{ 7 1+ ~ V-)- Moroever, at h l c = 1, m l c is a universal 
function of g c B /Tk; hence, the same is true for the in- 
frared singularity exponents r\ a . (At very large fields, 
however, a bulk Zeeman field, neglected above, will spoil 
universality, see Appendix 8.) 

For \B\ <C Tk, the magnetic moment is determined 
by the linear static susceptibility \o — 1/4Tk via m l c — 
-QeBxo, implying r\ a = rj - to lowest order in 

B/Tk- As \B\ increases past Tk and the magnetization 
tends towards m { c ~ — ^sgn(g c B), the exponents tend 
to r\ a — » — i(fig — l) 2 — <rsgn(g c B). Thus, they differ 
by 2 and attain opposite signs (since nf G [0,2]), de- 
pending on whether the photoexcited electron has spin 



a = =F S g n (.9c-B) (selectable by choice of circular polariza- 
tion of the incident ligth), i.e. whether its spin matches 
that of the lower or upper of the Zeeman-split e-levels, 
respectively. (Correspondingly, the photoexcited elec- 
tron will be said to have an "lower"- or "upper- level- 
spin", with exponent flower or Supper-) This implies a 
dramatic cr-dependence of the evolution of the lineshape 
A„{v) oc v^ 11 " with increasing \B\ (Fig. 5a). In the main 
text we focussed explicitly on the case that = 1, for 
which Slower/upper crosses over from | at B = to ±1 
for 1 2? | ^>T K (Eq. (9) and Fig. 5b): correspondingly the 
near-threshold singularity either becomes stronger, tend- 
ing towards v~ x ; or it becomes weaker, and once Supper 
turns negative, changes to an increasingly strong power- 
law decay, tending toward v +l . As mentioned in the 
main text, the difference between the two cases reflects 
the fact that in the limit \B\ ^> Tk, Anderson orthog- 
onality is absent or maximal for a = lower or upper, 
respectively. 

In this limit, the dependence of the lineshape on de- 
tuning and on a may also be interpreted in the following, 
alternative way, in terms of detuned transitions into the 
e-level: At the absorption edge, the photo-excitation of an 
upper-level-spin electron can be viewed as being a virtual 
excitation to the upper Zeeman level (with effective de- 
tuning given by ~ (v — g c B)) followed by a spin-flip of the 
e-level electron assisted by the creation of a FR electron- 
hole pair. Since the phase-space for creating these 
electron-hole pairs scales as v, the overall absorption rate 
for a = upper scales as vj (v — g c B) 2 ~ vj (g c B) 2 . In con- 
trast, for the photo-excitation of a lower-level-spin elec- 
tron, the effective detuning is is, and the excess energy 
is deposited into the Fermi sea via the creation of spin- 
preserving particle-hole pairs, for which the phase space 
again scales with v. This leads to an absorption rate for 
a = lower that scales as v/v 2 = \ jv. This explains the 
difference in power laws, flower — Supper = 2. 

The same picture works for larger photon energies 
above the edge, except that for the upper-level-spin elec- 
tron the energy deficit in the virtual state gets smaller. 
Once the energy deficit reaches zero, a resonant peak in 
the absorption results. Its width is determined by the 
Korringa relaxation rate for the upper e-level in the dot. 
In the scaling regime, the Korringa relaxation rate at 
large Zeeman energy (|p e -B| 3> Tk,T) is of the order of 
g e B/ ln 2 (g e J5/TK). Thus, the smaller the exchange inter- 
action with the band (i.e. the smaller Tk), the sharper 
the upper-level-spin absorption peak, and the smaller the 
weight outside the peak, near the absorption edge. 



APPENDIX 8: JUSTIFICATION OF 
NEGLECTING THE BULK MAGNETIC FIELD 

A magnetic field applied parallel to the Fermi reser- 
voir causes a Zeeman shift not only for the e- and h- 
levels, but also for the c-electrons in the Fermi reservoir, 
with g-factors of comparable magnitude: g c w —0.6-0.7, 
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FIG. S4: Dependence of magnetization on the bulk g- 
factor g c . The final magnetization is plotted as function 
of B/Tk for several different value of g c and two choices of 
g c : if g c = (dashed lines), the magnetization is sensitive 
to g c , but the crossover scale beyond which it approaches 
saturation is much larger than 7k. If g e = 2 (solid lines), 
we obtain the standard magnetization curve (also shown in 
Fig. 5b, dash-dotted line), for which the crossover scale is 7k. 
It's dependence on g c is so weak that it is not discernable 
unless magnified (see inset). NRG parameters: U = 0.17), 
e l e = 0.75E7, el = -0.517, V = 0.117, T K = 3.9 • 10" 2 I\ T = 0, 
A = 2.3, Kept states: 1024. 



<7h ~ 1.1-1.2 and g c « —0.2-0.4. Nevertheless, the effect 
of a bulk field, described by g c B^2 ka h(J(y kcr Cka, was ne- 
glected in the main text. The justification for this is as 
follows 31] : the effect of a bulk field on Kondo physics 
is to polarize the Fermi reservoir, causing the quantum 
dot to interact with a nonzero net spin (s£) w —pg c B/2. 
As a result, the effective Kondo coupling j^s c ■ s c gener- 
ates an extra contribution — (Jg c /2)Bs* to the local spin 
Hamiltonian, implying an effective local g-factor of 



.9, 



; ff = ,/, - g c J/2 . 



(S30) 



(We have checked that a similar result is obtained within 
the framework of the single-impurity Anderson model by 
using the scaling approach of Haldane to calculate the 
energy difference between having the local level occupied 
by a spin + or — .) Since the second term in Eq. (|S30[) in- 
volves the bare, not the renormalized exchange constant, 
its effect is small. 

A more detailed scaling analysis for \B\ ~^>Tk,v leads 
to the conclusion [3l| that for large magnetic fields the 
effective g-factor has the form 



i 



i 



2\n(B/T K ) 



The factor in square brackets describes a reduction of g° s 
with decreasing field that reflects the onset of the screen- 
ing correlations. In the so-called scaling limit where the 
coupling J is sent to zero and the bandwidth to infinity 



_ C3 
LU 




8000 1 0000 



B/T„ 



+ (,9c -,9c) J/2. (S31) 



FIG. S5: Dependence of the ground state energy on 
magnetic field. The final-state ground state energy (ground 
state energy of H i ), E f G , is plotted as a function of B /Tk as 
calculated by NRG (blue) and compared to the analytic ex- 
pressions given by Eq. (10) (red). Inset magnifies the region 
B <C Tk to show the quadratic dependence of the ground 
state energy (and hence the threshold) on the magnetic field. 
At large magnetic fields, the slope of the NRG curve dif- 
fers slightly from that of —\g e B due to slight (independent) 
renormalizations of the ecr level positions caused by its by 
hybridization with the FR. NRG parameters: U = 0.17), 
e* = 0.75(7, 4 = -0.5(7, T = 0.062(7, T K = 3.7- 10~ 3 Y, 
T — 0, A — 2.3, Kept states: 1024. 



in such a way that Tk remains fixed, the second term be- 
comes negligible for all fields smaller than the crossover 
field B* w Tke 1, /J ' . As a consequence, the local moment 
rn , calculated as a function of B, is expected to be in- 
dependent of g c for all \B\ < £>*; in other words, m e will 
be a universal function of g Q B in this regime. 

We have checked this expectation within the frame- 
work of a single-impurity Anderson model using NRG, 
proceeding as follows: The effect of a bulk field is to 
shift the Fermi reservoirs by —^(jg c B with respect to the 
Fermi energy at £p = 0. We thus used shifted band edges, 
^max/mm _ _ l^^g^ anc j em pi y ec i a cr-dependent 
Wilsonian logarithmic energy grid for each of the four 
energy intervals [0, D™ ax ] and [Df n ,0] (with a = ±). 
For the resulting Wilson chain we calculated the local 
moment m c (B) for both g c = and 2, for several differ- 
ent choices of g c , see Fig. EU As expected, the effect of 
turning on g c , though noticable for large fields if g c = 0, 
is very small if g c = 2. 
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